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Progress in lattice QCD 
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I review the recent progress in lattice QCD, which will be useful in heavy quark physics in the near future. 
Reviewing the theoretical developments in lattice QCD first, I focus our recent unquenched QCD with dynamical 
overlap fermion as implemented by JLQCD collaboration. I also introduce some of our recent studies on the 
B*Btt coupling and on the determination of \Vub\ through the dispersive bound. 
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1. Introduction 

The lattice computation of weak matrix elements 
can be defined as 



where 



lim 



a — ^O.rni — >m . 
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In order to extract physical matrix elements fi'om the 
lattice the sea quark effects, the renormalization fac- 
tor, and the continuum limit and chiral limit should 
be incorporated. In addition, the heavy quark should 
be treated either by the effective theory or by extrap- 
olation in the heavy quark mass from the smaller mass 
regime. These steps have been quite nontrivial tasks. 
Recently, there are three major progresses in lattice 
QCD which can drastically reduce the systematic er- 
rors in lattice QCD computation: (1) unquenched lat- 
tice QCD simulations in the chiral regime, (2) nonper- 
turbative renormalization, and (3) a new approach to 
heavy quarks on the lattice. In the following I explain 
these developments in some detail. 

1.1. Unquenched QCD simulations in the 
chiral regime 

A few years ago the only available large scale dy- 
namical QCD simulations with lightest pion mass 
~ 300 MeV, was by the MILC collaboration with 
the staggered quark [l| . This is because of the smallest 
numerical cost owing to the small degrees of freedom 
of the fermion and numerical stability from the ex- 
act partial chiral symmetry. In contrast, it has been 
thought that light dynamical fermion simulations in 
other fermion formalism would be difficult by Tfiops 
machines. 

In recent years, 0(10) Tfiops machines have become 
available in many places. Also the new precondition- 
ers for Dirac operator inversion algorithm such as 'do- 
main decomposition' [2;] and 'mass preconditioning' [3] 
enables us to treat the high- and low-mode contribu- 
tions to the Dirac operator separately. Combining 



this method with the 'multi-time scale' in the molec- 
ular dynamics step [1] gives significantly efficient al- 
gorithms for updating the gauge configurations in the 
hybrid Monte-Calro method. Owing to theses devel- 
opments, there are now may unquenched simulations 
|i, i, 0, B [i as shown in Table O Although the 
staggered fermion simulations is going much ahead, 
new results in other approaches will give important 
numerical and theoretical cross-checks with different 
advantages and disadvantages. 



1.2. Nonperturbative renormalization 

Another major development is the proposal of 
renormalization schemes with which lattice simula- 
tion can give the renormalization factors nonpertur- 
batively. 

One such scheme is the Schrodinger functional (SF- 
) scheme which is defined by amplitudes in a finite 
box with physical size L with Dirchlet boundary con- 
ditions in the temporal direction (lOj . The physical 
amplitudes to define this scheme are computable both 
perturbatively and nonperturbatively. The renormal- 
ization scale is defined as fj, — 1/ L. In this scheme 
both the renormalization constant of any local oper- 
ators and their running can be obtained. To obtain 
the running, additional simulations with different box 
sizes are needed. 

Another scheme is the regularization independent 
momentum scheme (RI-MOM) defined by off-shell 
quark/gluon amplitudes in Landau gauge The 
amplitudes are also computable both perturbatively 
and nonperturbatively. The renormalization scale is 
defined by the momentum scale. 



1.3. New approach to the heavy quark 

The precise computation of weak matrix elements of 
the B meson is one of the most important topics in lat- 
tice QCD. However since the typical lattice cutoff used 
in practical simulations is smaller than the bottom 
quark mass, naive lattice methods suffer from a large 
discretization error. For this reason lattice nonrela- 
tivistic QCD action has been widely used. Unfortu- 
nately, due to the nonrenormalizability of the action, 
one cannot take the continuum limit in this approach 
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Group 


Fermion Action 


Hf 


a(fm) 


L(fm) 


m^(GeV) 


MILC [IJ 


Improved Staggered 


2+1 


0.09, 0.12 


3 


> 300 


CERN [§.] 


Wilson, 0(a)-imp Wilson 


2 


0.052-0.075 


3 


> 300 


PACS-CS [6] 


0(a)-imp Wilson 


2+1 


0.07, 0.10, 0.12 


3 


> 210 


ETMC f7] 


twisted mass Wilson 


2 


0.075, 0.096 


3 


> 270 


RBC/UKQCD [8] 


Domain-wall 


2+1 


0.12 


3 


> 330 


JLQCD [9] 


Overlap 


2,2+1 


0.12 


2 


> 300 



Table I Projects of large scale lattice QCD simulations with light dynamical quarks. 



SO that one suffers from sizable systematic errors from 
discretization error and perturbative renormalization 
error. Recently new methods to treat the heavy-light 
meson in the continuum limit with nonperturbative 
accuracy have been proposed. 

The first approach was proposed by Alpha collab- 
oration [12] ■ They used the lattice HQET which is 
matched to QCD with nonperturbative accuracy by 
Schrodinger functional method in small volume with 
sufficiently fine lattice. Then, they evolve the lattice 
HQET action to coarser lattice by step scaling. They 
can also include l/Mh corrections into the action and 
operators. 

The second approach is the step scaling method pro- 
posed by the Rome II group [ij, The compute 
the physical observable (e.g. fs ) in small volume 
with Lq = 0.4 fm with sufficiently fine lattice using 
relativistic quark action. Then, they compute the fi- 
nite size corrections with larger volumes L = 2Lq, ALq 
by extrapolations from smaller heavy quark masses as 



(3) 



where (t(2Lo) 
Guazzini et al. 



/b,(2Lo) 



and a{Lao) 



[16j combined the above two meth- 
ods, i.e. they use the static limit result to interpolate 
the finite volume corrections. 

Table 11.31 shows the quenched QCD results of /s 
with nonperturbative accuracy. It is remarkable that 
all three approaches give consistent results with high 
accuracy. 



Method 


/b, (MeV) 


HQET with 1/M (Alpha) [15J 


193(7) 


Step scaling (Rome II) [14J 


195(11) 


Combination [161 


191(6) 



Table II Quenched QCD results of /s with nonperturba- 
tive accuracy. 



2. Recent results with dynamical overlap 
fermion 

The fermion action satisfying the Ginsparg- Wilson 
relation 



+ 75D = aD^oD (4) 
realizes the exact chiral symmetry on the lattice [3| 

5ip = 75(1 - aD)%l:, Stp = -ipj^. (5) 

An explicit construction of the Ginsparg- Wilson 
fermion called as the overlap fermion was proposed 
by Neuberger, which is defined as 



a 



(6) 



where Hw — l5{Dw + Mq) is the Wilson Hamiltonian 
with negative mass term Mq at the cutoff scale [T9| . 
The domain- wall fermion is another realization of the 
Ginsparg- Wilson fermion which introduce 5-th dimen- 
sion [20|. However, with finite extent in the 5-th di- 
mension, the chiral symmetry becomes approximate 
with an exponentially suppressed symmetry violation. 

The JLQCD collaboration succeeded in the first 
large scale lattice QCD simulation with a dynamical 
overlap fermion. Nf — 2 unquenched simulation on 
a 16^ X 32 lattice at a = 0.12 fm was carried out for 
6 quark masses covering the range m^/G rus and 
various measurements of physical observables were 

made [H, HE [H, M M M 

2.1. Chiral behavior of and /^r 

JLQCD collaboration studied the quark mass de- 
pendence of ttItt and /^r for Uf — 2 QCD with dynam- 
ical overlap fermion [9]. They found that the lattice 
data for < 450 MeV are wefi fitted with NLO 
ChPT formula 

rn^/niq = 2B[1 + x\nx) + cj,x (7) 
/tt = /(l-2a;lna;) + C4a;, (8) 

with X = ni^ / {Att f)'^ as shown in Figs[TJ They also 
studied the convergence of the ChPT by replacing the 
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expansion parameter x hy x = 2r7i^/(47r/)^ or ^ = 
ml / {At: f^f. They find that the NNLO ChPT with 
^-expansion can nicely describe the lattice data in the 
pion mass region of 290 ^ 750 MeV. 



2.2. B 



K 



Indirect CP violation in the K meson system ex 
is one of the most crucial quantities used to test the 
standard model and the physics beyond. The experi- 
mental value is determined with high accuracy as 



\tK\ = (2.233 ±0.015) X 10^3^ 



(9) 



Theoretically this quantity is described as |eA'| = 
f{p,r]) X C{fi) X Bk{p)- Here, f{p,r]) is a factor 
which depends on the CKM matrix elements, C(/i) is 
the Wilson coefficient from short-distance QCD cor- 
rections and Bk{p) is the bag parameter defined as 

R , , m - 75)^^7^(1 - 75)^] (A')!^"), ,^, 
Bk{p) = ^TJ-i (iO) 

S,JK"''K 

The main problem in unquenched lattice calculations 
is the possible operator mixing in wrong chiralities 
or tastes. The overlap fermion is free from opera- 
tor mixing owing to the exact chiral symmetry. The 
JLQCD collaborations study Bk with overlap fermion 
in 2 flavor QCD at lattice spacing a = 0.12 fm on 
physical volume with L = 2/m [l^]. They have 
4 points for the sea quark and 10 combinations of 
valence quark masses (toi, TO2)- They fit the data 
with NLO PQChPT. The renormalization factor is de- 
termined nonperturbatively by the RI-MOM scheme. 
They obtain B = 0.734(5)stat.(50)sys., where the dom- 
inant error comes from the finite size effect of order 
5%. 

It should be noted that the long standing operator 
mixing problem is solved with the advent of an overlap 
fermion with the exact chiral symmetry. Thus the 
above study heralds the beggining precision studies of 
Bk for which significant progress will be expected in 
near future. 



3. Some new results in heavy quark 
physics 

3.1. Determination of the B*Bti coupling 

The B*Bn coupling is a fundamental parameter of 
chiral effective Lagrangian with heavy-light mesons 
defined as 

L = -Tr [iiiv ■ DH] +.9bTr [HHA,, ■ 7^75] +0(1 /M), 

(11) 

where the low energy constant gj, is the B*B'k cou- 
pling, V is the four-velocity of the heavy-light meson 



are described by the B, 



r(l+7M«M)(^S75 + i?;7M), 



B or B*, and H, D 
B* and tt fields as H 
^ = exp(z^//), D^ = d^ + i(eta^^ + ^a^^t), = 

Once the B*Btt coupling is determined, the heavy 
meson effective theory can predict various quantities 
which are important for CKM phenomenology [3^ . 
For example the light quark mass dependence of the 
B meson decay constant can be determined as 



log(m2/A2 



(12) 

F is the low energy constant associated with the 
heavy-light axial- vector current. The form factor 
/+(g^) for the semileptonic decay B -kIv can also 
be expressed in terms of the B* meson decay constant 
fs' and cjb as 



2U 



9b 



k-A 



TUB 



(13) 



where v is the velocity of the B meson, k is the pion 
momentum, and A — tub* — tub- Therefore, the pre- 
cise determination of the B*B'k coupling is crucial for 
determining \ Vub\ and \ Vtd\ accurately. Despite its im- 
portance the B* Bn has been known not so accurately 
due to the large statistical error of the heavy-light me- 
son in static limit [H, [H, [11] . 

We carry out a precise determination of the B*BTr 
coupling in n/ = 2 QCD using 100 to 150 gauge con- 
figurations provided by CP-PACS collaboration (38| 
through JLDG (Japan Lattice DataGrid), which 12'^ x 
24 lattices at /3 = 1.80 and 16^ x 32 lattices at /3 = 1.95 
with two flavors of 0(a)-improved Wilson quarks and 
the Iwasaki gauge action. In order to improve the sta- 
tistical signal, we exploit the static quark action using 
the HYP smeared links 39] with the smearing param- 
eter values (ai,a2,a3) = (0.75^.6,0.3). We also em- 
ploy the all-to-all pro pag ator [40 ] using the low-mode 
averaging technique j4ll. 14^. Based on the previous 
quenched study [43], we take 200 low-eigenmodes. 

The physical value of the B*B'k coupling is obtained 
by multiplying the bare value by the renormaliza- 
tion constant at one-loop. The chiral extrapolation 
is made in three ways: (a) the linear extrapolation, 
(b) the quadratic extrapolation, and (c) the quadratic 
plus chiral log extrapolation where the log coefficient 
is determined from ChPT [4J]. We take the result at 
(3 = 1.95 as our best estimate for the physical value of 
cjoo, and estimate the discretization error of 0((aA)^) 
by order counting with A ~ 0.3 GeV. Including the 
perturbative error of 0{a'^) also by order counting, 
our result for is 



9o 



0.516(5)stat(31) chiral 

(28)pert(28) 



disc ■ 



(14) 
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Figure 1: ChPT fit witli x-expansion and xi-expansion. a;-expansion makes the convergence of the ChPT fits better. 
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Figure 2: The chiral extrapolation of the physical B* Bn 
coupling at j3 — 1.95. 



3.2. Dispersive bounds on the form 
factor for B nlu decay 

The momentum range of i? ^ ttIv form factors 
computed from lattice QCD is limited by the small 
recoil or large region. This leads to a big disad- 
vantage because most of the experimental data lies 
in large recoil region. While one can extrapolate in 

with a fit ansatz, this will always introduce some 
model dependence. Dispersive bounds is one possible 
way to constrain the q^ dependence in model indepen- 
dent fashion using unitarity. Consider the imaginary 
part of the vacuum polarization amplitude for the cur- 
rent V{x) — u'^^b{x) and a map as in Fig. 13.21 

! (15) 
Pl^"'{q) = i [ d'^xe"^ ''{0\T{V''{x)V'''<iO)}\0) 



Figure 3: A map from t plane to z plane 

= (9^9^ - g^-'q')n,{q^) + q^q'-Iloiq^m) 
Then, from dispersion relations one obtains 



1 92 



2a(q2)2 

d 



tlmlli{t) 



it 



with — —q^ and 77 an isospin factor, while x 's can 
be computed using the OPE and perturbative QCD. 
Unitarity tells us that this is equal to the sum over 
all the hadronic states, and dropping all the excited 
states and leaving only Bit and B* states gives an 
exact bound. 



[(t-t^)(t-t_)]3/2 

487r 



167r 



i3 



\F+{t)\^ < lmni(0, 
Fo{t)\^ < ImHoWClS) 



Combining Eqs. [T71[T51 and making change of variables 
in the integration from t to z 



z{t,to) 



(19) 
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with t± = {ms ± one obtains 

Wo|0/o)<Xo, (P<^/+|P0/+) <x+, (20) 

where J is a quantity which can be obtained using 
OPE and perturbative QCD. The inner product {g\h) 
for arbitrary functions g{z) and h{z) is defined by the 
integral along the unit circle in the z plane as 

{g\h)^ J ^Mz)r- (21) 

P{z) = z{t, rug) is multiplied by /+ in order to remove 
the B* pole inside the unit circle. Cauchy's theorem 
tells that if we know the additional integrated quantity 
{gi\P(j)+f+) with a set of known functions {gi{z),i = 
1, iV} one can make the bound stronger as 



found that the chiral behavior of to^ and are consis- 
tent with Next-to-next-leading order Chiral Perturba- 
tion theory. With the advent of the exact chiral sym- 
metry a precise determination of BK was discussed. 

I also explained some new results in heavy quark 
physics such as B* Bi: coupling and the model inde- 
pendent determination of \Vub\ from the dispersive 
bound. 
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stronger form factor bounds with statistical analy- 
sis. We improved the bound using also the experi- 
mental spectrum from CLEG as additional inputs 
[4^ . After the B AB AR measurement of the spec- 
trum of -B ^ ttIv decay [Jt'I, Arnesen et al. [4^ 
set gn{z) = z" to obtain a simple bound on the co- 
efficients of the z-polynomial parameterization, which 
was further improved imposing HQET power counting 
by Becher and Hill [49], [50]. 

Since the BABAR measurement of the spectrum 
allows for the form factor shape determination, we 
also updated our determination of the \Vub\ using the 
dispersive bound [5l[. Using the form factor from 
HPQCD collaboration [13] and the CLEG data, we 
obtain our preliminary estimate 



\Vub 



[3.41°;^] X 10 



-3 
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4. Summary 

There has been major progress in the unquenched 
QCD simulation in the chiral regime, the renormal- 
ization schemes which allows for nonperturbative de- 
terminations of the renormalization factors, and in 
the new approach to the heavy quarks on the lattice. 
These developments have been tested in light hadron 
physics or in quenched QCD and are promising for 
improving the lattice calculation for B physics in the 
near future. 
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overlap fermion by the JLQCD collaboration. It was 
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